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^ Abstract 

I review recent works on the symmetry and the structure of the superconducting gap in Fe-based 

(N 

^ superconductors and on the underlying pairing mechanism in these systems. The experimental data 

O 

on superconductivity show very rich behavior, with potentially different symmetry of a supercon- 

ducting state for different compositions of the same material. The variety of different pairing states 

f3 raised the issue whether the physics of Fe-based superconductors is model-dependent or is univer- 
O 

O sal, governed by a single underlying pairing mechanism. I argue that the physics is universal and 
U 

Qh that all pairing states obtained so far can be understood within the same universal pairing scenario 

^ and are well described by the effective low-energy model with a small number of input parameters. 
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I. INTRODUCTION 



The discovery, in 2008, of superconductivity in Fe-based pnictides [1] (binary compounds 
of the elements from the 5th group: N, P, As, Sb, Bi) was, arguably, among the most 
significant breakthroughs in condensed matter physics during the past decade. A lot of 
efforts by the condensed-matter community have been devoted in the few years after the 
discovery to understand normal state properties of these materials, the pairing mechanism, 
and the symmetry and the structure of the pairing gap. 

The family of Fe-based superconductors (FeSCs) is already quite large and keeps growing. 
It includes various Fe-pnictides such as 1111 systems RFeAsO {R —rare earth element) [1-4], 
122 systems XFe2As2(X=alkaline earth metals) [5-7], 111 systems hke LiFeAs [9], and also 
Fe-chalcogenides (Fe-based compounds with elements from the 16th group: S, Se, Te) such 
as FeTei„^Sc^ [10] and AFe.Ses (A = K, Rb, Cs) [11,12]. 

Parent compounds of FeSCs are metals, in distinction to cuprate superconductors for 
which parent compounds are Mott insulators. Still, in similarity with the cuprates, in most 
cases these parent compounds are antiferromagnetically ordered [13]. Because electrons 
which carry magnetic moments still travel relatively freely from site to site, the magnetic 
order is often termed as a spin-density- wave (SDW), by analogy with e.g., antiferromagnetic 
Cr, rather than "Heisenberg antiferromagnetism" - the latter term is reserved for systems 
in which electrons are "nailed down" to particular lattice sites by very strong Coulomb 
repulsion. 

Superconductivity (SC) in FeSCs emerges upon either hole or electron doping (see Fig. 
1), but can also be induced by pressure or by isovalent replacement of one pnictide element 
by another, e.g.. As by P (Ref. [8]). In some systems, like LiFeAs [9] and LaFePO [14], SC 
emerges already at zero doping, instead of of a magnetic order. 

The magnetism, the electronic structure, the normal state properties of FeSCs, and the 
interplay between FeSCs and cuprate superconductors have been reviewed in several recent 
publications [15-25]. This review is an attempt to summarize our current understanding 
of the pairing mechanism and the symmetry and the structure of the pairing gap at various 
hole and electron dopings. 

The phenomenon of SC has a long history. SC has been discovered by Kamerlingh Onnes 
exactly a century ago [26]). It has been explained in general terms nearly fifty years later, in 
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FIG. 1: Schematic phase diagram of Fe-based pnictides upon hole or electron doping. In the shaded 
region, superconductivity and antiferromagnetism co-exist. Not all details/phases are shown. Su- 
perconductivity can be initiated not only by doping but also by pressure and/or isovalent replace- 
ment of one pnictide element by another [8]. Nematic phase at T > is subject of debates. 
Superconductors at large doping are KFe2As2 for hole doping [40,41] and A^Fe2-ySe2 (A = K, Rb, 
Cs) for electron doping [11,12]. Whether superconductivity in pnictides exists at all intermediate 
dopings is not clear yet. From Ref. [24]. 

1957, by Bardeen, Cooper, and Schrieffer (BCS), who demonstrated that an arbitrary weak 
attractive interaction between low-energy fermions is sufficient to pair them into a bound 
state. At weak coupling, paired fermions immediately form Bose-Einstein condensate and 
behave as one single macroscopic quantum object and move coherently under the applied 
electric field, i.e superconduct. In (i-dimensional electronic systems low-energy fermonic 
states are located, in momentum space, near particular d — 1 dimensional surfaces, called 
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Fermi surfaces (FS) on which fermionic energy is zero relative to the chemical potential. 
At weak/moderate coupling, the pairing problem is confined to a near vicinity of a FS. 
The interaction between fermions is generally non-singular with respect to variations of the 
distance to the FS and can be approximated by its value right on the FS. 

What causes the attraction between fermions is a more subtle question, and the nature 
and the origin of the pairing glue have been the subject of great debates in condensed-matter 
community over the last 50 years. BCS attributed the attraction between fermions to the 
underlying interaction between electrons and phonons [27] (the two electrons effectively 
interact with each other by emitting and absorbing the same phonon which then serves as 
a glue which binds electrons into pairs). Elcctron-phonon mechanism has been successfully 
applied to explain SC in a large variety of materials, from Hg and Al to recently discovered 
and extensively studied MgB2 with the transition temperature Tc — 39K [28]. Non- 
phononic mechanisms of the pairing have also been discussed, most notably in connection 
with superfiuidity in ^He [29] , but didn't become the mainstream before the discovery of SC 
in LaBaCuO in 1986 [30]. That discovery, and subsequent discoveries of superconductivity 
at higher Tc in other cuprates signaled the beginning of the new era of "high-temperature 
superconductivity" to which FeSCs added a new avenue with quite high traffic over the last 
three years. 

Superconductivity is quite robust phenomenon. It has been known from early 60th [31] 
that in isotropic systems the equation for superconducting factorizes if one expands the 
interaction between the two fermions in partial components corresponding to interactions in 
the subspaces with a given angular momentum of the two interacting fermions I = 0, 1, 2, 3, 
etc [in spatially isotropic systems / = component is called s— wave, / = 1 component is 
called p— wave, I — 2 component is called d— wave, and so on]. If just one component with 
some I is attractive, the system undergoes a SC transition at some temperature T — T^.. For 
phonon-mediate superconductors, s— wave superconductivity is the most likely outcome. In 
the cuprates, however, the pairing symmetry has been firmly established as d— wave. The 
vast majority of researches believe that such pairing is not caused by phonons and emerges 
instead due to screened Coulomb interaction between electrons. The screened Coulomb 
interaction U (r) is constant and repulsive at short distances but has a complex dependence 
on r at large distances and may develop an attractive component at some I. One solid reason 
for the attraction, at least at large I, has been identified by Kohn and Luttinger back in 
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1965 (Ref. [32]). 

In lattice systems, angular momentum is no longer a good quantum number, and the 
equation for T^. only factorizes between different irreducible representations of the lattice 
space group. In tetragonal systems, which include both cuprates and FeSCs , there are four 
one-dimensional irreducible representations Aig, Big, B2g, and and one two-dimensional 
representation E2g. Each representation has infinite set of eigenf unctions. The eigen- 
functions from Aig are invariant under symmetry transformations in a tetragonal lattice: 
X — > —X, y —y, X — >■ y, the eigenfunctions from Big change sign under x — >• y, and so 
on. If a superconducting gap has Aig symmetry, it is often called s— wave because the first 
eigenfunction from Aig group is just a constant in momentum space (a 5— function in real 
space). If the gap has Big or B2g symmetry, it is called li— wave {d.j,2_y2 or dxy, respectably), 
because in momentum space the leading eigenfunctions in Big and B2g are cos kx — cos ky 
and sin k^ sin ky, respectively, and these two reduce to I — 2 eigenfunctions cos 26 and sin 26 
in the isotropic limit. 

In the cuprates, the superconducting gap has been proved experimentally to have Big 
symmetry [33]. This gap symmetry appears quite naturally in the cuprates, in the doping 
range where they are metals, if one assumes that the glue that binds fermions together is 
a spin- fluctuation exchange rather than a phonon (see Fig. 2). The notion of a spin fluctu- 
ation is actually nothing but the convenient way to describe multiple Coulomb interactions 
between fermions. It is believed, although not proved rigorously, that in systems located 
reasonably close to a magnetic instability, the fully screened Coulomb interaction between 
fermions can be approximated by an effective interaction in which fermions exchange quanta 
of their collective fluctuations in the spin channel. That Big gap is selected is not a surprise 
because such gap A(A;) oc coskx — cos ky changes sign not only under k^ ky but also 
between k and k' = k -|- Q where Q = (tt, tt) is the momenta at which spin fluctuation- 
mediated pairing interaction C/(k, k') is peaked. This sign change is the crucial element for 
any electronic mechanism of superconductivity because one needs to extract an attractive 
(negative) component from repulsive (positive) screened Coulomb interaction. For Big gap 
such a component is / (ik(ik'A(A;)f/(k, k')A(A;'), and the integral obviously has a negative 
value when t/(k, k') is peaked at (tt, tt). 

In FeSCs, magnetism and superconductivity are also close neighbors on the phase dia- 
gram, and it has been proposed [34,35] at the very beginning of the Fe era that the pairing 
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FIG. 2: A comparison of the pairing state from spin-fluctuation exchange in cuprate SCs and in 
FeSCs. In the cuprates (left panel) the FS is large, and antiferromagnetic Q = (vr, vr) connects 
points on the same FS. Because spin-mediated interaction is positive (repulsive), the gap must 
change sign between FS points separated by Q. As the consequences, the gap changes sign twice 
along the FS. This implies a d— wave gap symmetry. In FeSCs (left panel) scattering by Q moves 
fermions from one FS to the other. In this situation, the gap must change sign between different 
FS, but to first approximation remains a constant on a given FS. By symmetry, such a gap is an 
s— wave gap. It is called s'' because it changes sign between different FSs 

mechanism in FeSCs is also a spin-fluctuation exchange. However, the geometry of low- 
energy states in FeSCs and in the cuprates is different, and in most FeSCs the momentum 
Q connects low-energy fermionic states near the center and the corner of the Briilouin zone 
(see Fig.2). A simple experimentation with trigonometry then tell us that the SC gap A{k) 
must be symmetric with respect to ky and kx —kx, but still must change sign under 
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k — >■ k + Q. Such gap belongs to Aig representation, but it only has contributions from a 
particular subset of Aig states with the form cosk^ + cos ky, cos3A;^ + cos{3ky), etc which 
all change sign under k — >■ k + Q. Such gap is generally called an extended s— wave gap, or 
s'^" gap. 

Majority of researches do believe that in weakly/moderately doped FeSCs the gap does 
have s'^" symmetry However, numerous studies of superconductivity in FeSCs over the last 
three years demonstrated that the physics of the pairing is more involved than it was origi- 
nally thought because of multi-orbital/multi-band nature of low-energy fermionic excitations 
in FeSCs (see below) . It turns out that both the symmetry and the structure of the pairing 
gap result from rather non-trivial interplay between spin-fluctuation exchange, intraband 
Coulomb repulsion, and momentum structure of the interactions. In particular, an s^wave 
gap can be with or without nodes, depending on the orbital content of low-energy excita- 
tions. In addition, the structure of low-energy spin fluctuations evolves with doping, and 
the same spin-fluctuation mechanism that gives rise to s'^~ gap at small/moderate doping 
in a particular material can give rise to a d— wave gap at strong hole or electron doping. 

There is more uncertainly on the theory side. In addition to spin fluctuations, FeSCs also 
possess charge fluctuations whose strength is the subject of debates. There arc proposals [36, 
37] that in multi-orbital FeSCs charge fluctuations are strongly enhanced because the system 
is reasonably close to a transition into a state with an orbital order (e.g., a spontaneous 
symmetry breaking between the occupation of different orbitals). A counter-argument is 
that orbital order does not develop on its own but is induced by a magnetic order [38]. If 
charge fluctuations are relevant, one should consider, in addition to spin-mediated pairing 
interaction, also the pairing interaction mediated by charge fluctuations. The last interaction 
can give rise to a conventional, sign-preserving s— wave pairing [36]. A "p-wave" gap scenario 
(a gap belonging to representation) has also been put forward [39] . 

From experimental side, s-wave gap symmetry is consistent with ARPES data on mod- 
erately doped KFe2As2 and BaFe2(Asi_a;Pa;)2, which detected only a small variation of the 
gap along the FSs centered at (0,0) (Ref. [42]), and with the evolution of the tunnehng 
data in a magnetic fleld [98. However, for heavily hole-doped KFc2As2various experimental 
probes [41] indicate the presence of gap nodes, which for the FS geometry in these materi- 
als [40] are consistent with a d-wave gap. For the doping range where the gap is very likely 
an s— wave, the data on some FeSCs were interpreted as evidence for the full gap [44-47, 
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while the data for other FeSCs were interpreted as evidence that the gap has nodes [48,49] or 
deep minima [50-52]. In addition, recent nuclear magnetic resonance (NMR) experiments 
on LiFeAs have been interpreted in favor of a p-wave gap [43] . 

In this paper, I argue that all these seemingly very different gap structures (with the 
exception of a p— wave), actually follow quite naturally from the same underlying physics 
idea that FeSCs can be treated as moderately interacting itinerant fermionic systems with 
multiple FS sheets and effective four-fermion intra-band and inter-band interactions in the 
band basis. I introduce the effective low-energy model with small numbers of input param- 
eters [66] and use it to study the doping evolution of the pairing in hole and electron-doped 
FeSCs. It has been argued [66] that various approaches based on underlying microscopic 
model in the orbital basis reduce to this model at low energies. 

The paper is organized as follows. In Sec. II I discuss general aspects of the pairing in 
FeSCs. I briefly review the band structure of FeSCs and show that it contains several bands 
of low-energy excitations. I then present generic symmetry considerations of the pairing in 
a multi-band superconductor. I show that a "conventional wisdom" that an s-wave gap is 
nodeless along the FSs, d-wave gap has 4 nodes, etc, has only limited applicability in multi- 
band superconductors, and there arc cases when the gap with four nodes has an s— wave 
symmetry, and the gap without nodes has a o?— wave symmetry. In Sec. Ill I discuss the 
interplay between intra-band and inter-band interactions, first for a toy two-pocket model 
and then for realistic multi-pocket models, and set the conditions for an attraction in an 
s— wave or a d— wave channel. I consider 5-orbital model with local interactions, convert 
it into a band basis, and argue that for most of input parameters the bare interaction is 
repulsive in all channels due to strong intra-pocket Coulomb repulsion. In Sec. IV I discuss 
the ways to overcome Coulomb repulsion. I review random phase approximation (RPA) and 
renormalization group (RG) approaches and show that magnetic fluctuations enhance inter- 
pocket interaction, if this interaction is positive, what gives rise to an attraction in both 

and dx^^y2 channels. I briefly discuss s"*""^ pairing, which emerges if input parameters 
are such that inter-pocket interaction is negative. In Sec.V I use the combination of RPA 
and leading angular harmonic approximation (LAHA) to analyze the pairing in s— and 
(i— wave channels at different dopings. I show that magnetically-mediated pairing leads to 
(i) an s"^ superconductivity with nodes on electron pockets for moderate electron doping, (ii) 
an 5=*= superconductivity without nodes for moderate hole dopings, (iii) a nodeless d— wave 
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FIG. 3: The electronic structure of FeSCs. In weakly and moderately electron-doped materials 
(left panel) the FS consists of quasi-2D warped cylinders centered at (0, 0) and (vr, vr) in a 2D 
cross-section. The ones near (0,0) are hole pockets (filled states are outside cylinders), the ones 
near (vr, vr) are electron pockets (filled states are inside cylinders) There also exists a quasi-3D hole 
pocket near = vr. In hole-doped FeSCs the electronic structure is very similar, but 3D hole 
pocket becomes quasi-2D warped hole cylinder. From Ref. [23]. 

superconductivity for strong electron doping (except, possibly, small nodal regions near 
kz = and (iv) a nodal d— wave superconductivity for strong hole doping. I briefly 

review the experimental situation in Sec. VII and present concluding remarks in Sec. VII. 
I list separately the summary points and the list of future issues. 

II. GENERIC ASPECTS OF PAIRING IN FESCS 
A. The electronic structure 

The electronic structure of FeSCs at low energies is rather well established by ARPES [53] 
and quantum oscillation measurements [54]. In weakly and moderately electron-doped ma- 



9 



terials, like BaFei_a;Co^Fe2As2 the FS contains several quasi-2D warped cylinders centered 
ai k — (0, 0) and k — (tt, tt) in a 2D cross-section, and may also contain a quasi-3D pocket 
near k^ — n (Fig.3). The fermionic dispersion is electron-like near the FSs at (7r,pi) (filled 
states are inside a FS) and hole-like near the FSs centered at (0, 0) (filled states are outside 
a FS). In heavily electron-doped FcSCs, like Aa;Fei_ySe2 (A = K, Rb, Cs), only electron 
pockets remain, according to recent ARPES studies. [11] In weakly and moderately hole- 
doped FeSCs, hke Bai_^K3;Fe2As2, the electronic structure is similar to that at moderate 
electron doping, however the spherical FS becomes the third quasi 2D hole FS centered at 
(27r, 0) = (0,0). In addition, new low-energy hole states likely appear around {n,n) and 
squeeze electron pockets [55]. At strong hole doping, electron FSs disappear and only 
only hole FSs are present [40] These electronic structures agree well with first-principle 
calculations [17,56,57], which is another argument to treat FeSCs as itinerant fermionic 
systems. The measured FS reflects the actual crystal structure of FeSCs in which there are 
two non-equivalent positions of a pnictide above and below an Fe plane, and, as a result, 
there are two Fe atoms in the unit cell (this actual Brillouin zone (BZ) is called "folded 
BZ"). From theory perspective, it would be easier to work in the BZ which contains only 
one Fe atom in the unit cell (this theoretical BZ is called "unfolded BZ"). I illustrate the 
difference between folded and unfolded BZ in Fig.4. In general, only folded BZ is physically 
meaningful. However, if by some reason a potential from a pnictogen (or chalcogen) can be 
neglected, the difference between the folded and the unfolded BZ becomes purely geomet- 
rical: the momenta k^ and ky in the folded BZ are linear combinations of k^ and ky in the 
unfolded BZ: k^ = kx + ky, ky = kx — ky. In this situation, folded and unfolded BZ become 
essentially equivalent. 

Most of the existing theory works on the pairing mechanism and the structure of the 
SC gap analyze the pairing problem in the unfolded BZ, where which two hole pockets are 
centered at (0,0) and one at (7r,pi), and the two electron pockets are are at (0, tt) and 
(tt, 0). It became increasingly clear recently that the interaction via a pnictogen/chalcogen 
and also 3D effects do play some role for the pairing, particularly in strongly electron-doped 
systems. [58,59] However, it is still very likely that the key aspects of the pairing in FeSCs 
can be understood by analyzing a pure 2D electronic structure with only Fe states involved. 
Below I assume that this is the case and consider a 2D model in the unfolded BZ with hole 
FSs near (0,0) and {n,pi) and electron FSs at (0, tt) and (7r,0). 
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FIG. 4: Upper panel: 3D electronic structure of LaOFeAs (left) and its 2D cross-section (left). In 
only Fe states are considered, an elementary cell contains one Fe atom (green). The actual unit 
cell (blue) contains two Fe atoms because of two non-equivalent positions of a pnictide above and 
below the Fe plane. Lower panel - the location of hole and electron FSs in a 2D cross section in 
the folded BZ (two Fe/cell, right) and in the unfolded BZ (one Fe/cell, left). From Refs. [60], 
[61] (b) and [62]b. 

B. The structure of s— wave and d— wave gaps in a multi-band SC 

I now use the multi-band electronic structure as input and consider the pairing problem 
at weak coupling. I show that an s— wave gap generally has angle dependence on electron 
FSs and may even has nodes, while a d-wave gap which is normally assumed to have nodes, 
may in fact be nodeless on the same electron FSs. 



11 




FIG. 5: Intra-pocket and inter-pocket interactions in 4-band 2D model for FeSCs. For simplicity, 
only one hole FS is shown. F, X, and Y points are (0, 0), (tt, 0), and (0,7r), respectively. is 
the interaction within hole pocket, Fej/i and Fgjft are the interactions between hole and one of 
electron pockets, and Fejen ^e2,e2 (not shown) and Gammaeie2 are intra-pocket and inter-pocket 
interactions involving the two electron pockets. Each interaction contains s— wave and d— wave 
component, and even for s— wave T^.h and Fe-e^ depend on the angles along electron FSs. From 
Ref. [64]. 

A generic low-energy BCS-type model in the band basis is described by 

'H = Y^ ei(k)aL«ik + ^ijlk, p)a!k«I-k«ipaj-P (1) 

j,k ij,k,p 

The quadratic term describes low-energy excitations near hole and electron FSs, labeled by 
i and j, and the four-fermion term describes the scattering of a pair {k f, —k l) on the FS 
i to a pair {p f, —p l) on the FS j. These interactions are either intra-pocket interactions 
(hole-hole F/^-zj. or electron-electron Fg^ei ), or inter-pocket interactions (hole-electron T^jhi, 
hole-hole Thij^hj, and electron-electron Te^^ej)- I illustrate this in Fig.5. 

Assume for simplicity that the frequency dependence of F can be neglected and low- 
energy fermions are Fermi-liquid quasiparticles with Fermi velocity Vkp- In this situation, 
the gap A(fc) also doesn't depend on frequency, and the linearized gap equation becomes 
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the eigenf unction/eigenvalue problem: 

A,A,(^) = - 1 |^r(k^, p^) A,(p) (2) 

where Aj are eigenfunctions and Aj are eigenvalues. The system is unstable towards pairing 
if one or more Aj are positive. The corresponding Tc^j scale as Tc^i — AjC"^/'*'*. Although Aj 
are generally different for different i, the exponential dependence on 1/Aj implies that, most 
likely, the solution with the largest positive Aj emerges first and establish the pairing state, 
at least immediately below Tc. 

Like wc discussed in the Introduction, the vertex r[k,p) can be decomposed into rep- 
resentations of the tetragonal space group (one-dimensional representations arc Aig, Big, 
B2g, and ^2g)- Basis functions from different representations do not mix, but each contains 
infinite number of components. For example, s— wave pairing corresponds to fully symmet- 
ric Aig representation, and the s— wave (Aig) component of r{k,p) can be quite generally 
expressed as 

r(^^)(fc,p) = TsiKp) = E^mn*^(^)*n(p) (3) 

m,n 

where ^^(/c) are the basis functions of the Ai^ symmetry group: 1, cos 

kx cos ky ^ cos kx + 

cosky, etc, and A^^^ are coefficients. Suppose that k belongs to a hole FS and is close to 
= 0. Expanding any wave function with Aig symmetry near A; = 0, one obtains along 
|k| = kp, 

*^(A;) = + bmCos44>k + c^cosS^^ + ... (4) 

where (pk is the angle along the hole FS (which is not necessary a circle) . Similarly, for Big 
representation the wavefunctions are coskx — cos ky, cos2kx — cos2ky, etc, and expanding 
them near k = one obtains 

^m(^) = «m COS 2(pk + b*^ COS 60fe + C*^ COS lO0fc + ... (5) 

There are no fundamental reasons to expect that bm, Cm or b^,d^ are much smaller than 
Gm or aj^, but sub-leading terms are often small numerically. The known example is the 
numerical smallness of cos(2 + 4n)0 components with n > 1 of the (ij.2_^2— wave gap in the 
spin-fluctuation scenario for the cuprates [63] (meaning that a d-wave gap can be reasonably 
well approximated by cos 20). Taking this example as circumstantial evidence, I neglect all 
subleading terms, i.e., assume that s— wave interaction between fermions on the hole FSs can 



13 



be approximated by an angle-independent ri.f^.{k,p) — Uhfhj {hi label different hole FSs), 
while d— wave (Big) interaction can be approximated by r^.^^(A;,p) = Uhihj cos20fecos20p. 

The situation changes, however, when we consider the pairing component involving 
fermions from electron FSs. Suppose that k are still near the center of the BZ, but p 
are near one of the two electron FSs, say the one centered at (0,7r). Consider all possible 
^n(p) with Aig symmetry A simple experimentation with trigonometry shows that there are 
two different subsets of basis functions: 

subset 1:1, cos cos Py, cos2pa. + cos2pj,... 

subset II : cospx + cospy, cos Spx + cos Spy... (6) 

For a circular FS centered at (0,7r), the functions from subset 1 can be again expanded in 
series of cos4Z0p with integer /. The functions from subset II are different - they all vanish 
at (0, tt) and are expanded in series of cos(20p + 4Z0p) (the first term is cos 20^, the second 
is cos60p, etc). For elliptic FS cos4Z0p and cos(20p + 4Z0p) terms appear in both subsets. 
In both cases, the total 

*m(p) = am + bm COS 40j, + Cm COS 8(j)p + ... 

+am COS 20p + bm cos 60p + cos 100^ + ... (7) 

For the other electron FS, ^^(p) is the same, but momentum components p^ and Pj^are 
interchanged, hence the sign of all cos(20 + 4Z0p) components changes. 

Let's make the same approximation as before and neglect all components with I > 0. 
Then 

*m(p) = Qm ± Om COS 2(pp (8) 

where the upper sign is for one electron FS and the lower for the other. It is essential that 
the angle-independent term and the cos 20p term have to be treated on equal footing because 
each is the leading term in the corresponding series. Combing (8) with the fact that "iffnik) 
can be approximated by a constant, we obtain a generic form of the s— wave component of 
the interaction between fermions near hole and electron FSs 

KuhSf'^P) = t/e,,,,(l + 2ae,hcos20p^, + ...) 

K,x(k,P) = Uex (1 - 2ae,/.cos20p^, + ....) (9) 
where dots stand for cos40fe, cos4(/)p, cos60p, etc terms. 
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By the same reasoning, s— wave components of inter-pocket and intra-pocket interactions 
between fermions from electron FSs are 

Kueii^^P) = Ue,e (1 + S^ee (cOS 20fe^, + COS 20^^ J 

+4/3ee COS 2(j)k,i COS 2(pp^^ + ... 
K2,e2i^^P) = Ue,e (1 " ^^ee (cOS 20fc^2 + COS 20^^ J 

+4:/3ee COS 20fc^2 COS 20p^2 + ■■■ 
Ki,e2i^^P) = ^e,e (1 + 2aee (Cos20fe^l - COs20p^2) 

-4/3ee COS 20fc^j COS 20p^2 + ... (10) 

Once the pairing interaction has the form of Eqs. (9) and (10), the gaps along the hole 
FSs are angle-independent (modulo cos 40 terms) , but the gaps along the two electron FSs 
are of the form 

A(^)(A;) = Ae±AeCos20fc. (11) 

When Ae is small compared to Ag, the angle dependence is weak, but when |Ae| > |Ae|, 
s— wave gaps have nodes at "accidental" values of 4>, which differ between the two electron 
FSs. 

A similar consideration holds for (1^2 _y2 gap. Within the same approximation of leading 
angular momentum harmonics, we have 

rei,ft,(^,P) = C/e,/,,COs20ft. (l + Q;e,feCOS20p^J + ... 

re2A(^,P) = C/e,/.,cos20^,(-l + ae,/.cos20p^J + ... (12) 

and 

rl,ei(^>P) = Ue,e (1 + 2aee (cOS 20^^, + COS 20^^ J 

+4:/3ee COS 20fcei COS 2(f)p^^ + ... 
1^62,62 (^'^^) = ^e,e (1 - 2aee {cOs2(j)k^^ + COS 2 J 

+4:/3ee COS 20jfc^2 COS 2(f)p^^ + ... 
K,e2ik,p) = Ue,e ("1 - 2o;ee (cOS 20fe^^ - COS 2 J 

+4:/3ee COS 24>k^-^ COS 20^^^ + •■■ (13) 
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The solution of the gap equation then yields the gap in the form 

A'^f'^ (k) = ±Ae + Ae COS 20fe. (14) 

Along the hole FS, the gap behaves as a conventional d— wave gap with 4 nodes along the 
diagonals. Along electron FSs, the two gaps differ in the sign of the angle-independent 
terms, and have in- phase cos 20 oscillating components. When Ag « A^ the two electron 
gaps are simply "plus" and "minus" gaps, but when Ag > Ag, each has accidental nodes, 
again along different directions on the two electron FSs. 

We see therefore that the geometry of the FSs in FeSCs affects the gap structure in 
quite fundamental way: because electron FSs are centered at the k points which are not 
along BZ diagonals, s— wave gaps on these FSs have cos 2(f) oscillations which one normally 
would associate with a d— wave symmetry, and d— wave gaps have constant (plus-minus) 
components which one would normally associate with an s— wave symmetry. When these 
"wrong" components are large, the gaps have accidental nodes. These nodes may be present 
or absent for both s— wave and d— wave gaps, i.e., symmetry constraints play no role here. 

An s— wave gap with nodes in one of the "exotic" options offered by the electronic struc- 
ture of FeSCs. Another "exotic" option is a d— wave state without nodes. In heavily electron- 
doped FeSCs, hole states are gapped, and only electron FSs remain. The d— wave gaps on 
these two FSs have no nodes if cos 2(f) oscillation component is smaller than a constant term, 
hence the system will display a behavior typical for a fully gapped SC despite that the gap 
actually has a d-wave symmetry. There are even more exotic options offered by the actual 
three-dimensionahty of the electronic structure and/or the hybridization of the electron FSs 
due to interaction via a pnictide/chalcogen. One RPA calculation [59] places the nodes of 
the gap on hole pockets, near particular k^- Another study argues [58] that a d— wave 
gap for heavily electron doped FeSCs must have nodes near = re/ 2. 

It is essential, however, that either a more conventional one or a more exotic pairing 
state develops only if the corresponding eigenvalue is positive. To understand under what 
conditions Aj > we now have to consider specific models for FeSCs. 
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III. THE INTERPLAY BETWEEN INTRA-POCKET AND INTER-POCKET IN- 
TERACTIONS 

In this section we make the first step in the analysis of what causes the attraction in 
FeSCs and consider how the sign and magnitude of Aj depends on the interplay between 
intra-pocket and inter-pocket interactions. 



A. Toy two-pocket model 

As a warm-up, consider first an idealized two-pocket model consisting of two identical 
circular pockets: a hole pocket at (0, 0) and an electron pocket at (tt, tt). Because the electron 
FS is centered at k along the diagonal, the cos 2(f) terms in the electronic gap are no longer 
present, so some part of the physics of FeSCs is lost. Still, this idealized model is a good 
staring point to consider the interplay between intra-pocket and inter-pocket interactions. 

Compared to Eqs. (9-13), we now have T^f^ = Uhh, = Uhe-i F^^ = t/ge, and 

^hhikp^PF) = Uhh cos 2(l)k cos 2(j)p, ^heik.p) = C/fteCOs26'fe, F^g = C/ee. 

The eigenvalue problem, Eq. (2), reduces to the set of two coupled equations in each 
channel. Solving them, we obtain 



_ -{Uhh + Uee) ± \J {Uhh - UeeY + 4^1^ 

Ai,2 — ^ 



and 



^^1,2 — V-^) 



where — UijNp and Np — m/(27r) is the density of states (DoS) per spin projection. 

We see that a solution with a positive A in either s— or d— wave channel exists if the 
inter-pocket interaction is larger than intra-pocket interactions. Specifically, one needs 

> UeeUhh or > UeeUhh- (17) 

If this condition is not met, the system remains in the normal state down to T = 0. If 
Eq. (17) is satisfied and Uhe > 0, s— wave solution with A* > yields a "plus-minus" gap, 
Ag = — A/j («=*= state). If Uhe < (what requires intra- band attraction), s— wave solution 
is a conventional one, with Ag = Ah- Similarly, d— wave solutions yield cos 20 gaps on the 
hole and electron FSs with either zero or tt phase shift. 

17 



B. Multi-band models 



Consider next a more realistic case of two electron FSs centered at (0, tt) and (tt, 0). 
Now hole-electron and electron-electron interactions have cos 20 terms, and the eigen- 
value/eigenf unction problem, Eq. (2), reduces to the set of either four (or five) coupled 
equations in either s-wave or d-wave channels: two (or three) A's are the gaps on the hole 
FSs, and two other A's are angle-independent and cos 20 components of the gaps on the 
electron FSs. Accordingly, there are either four or five different and A,^. 

The analysis of 4 x 4 or 5 x 5 gap equations is tedious but straightforward. I will not 
discuss it in length (for a detailed discussion see Refs. [64,66]) but rather focus on an issue 
of whether it is still required that the inter-pocket interaction Uhe exceeds the threshold 
set by intra-pocket hole-hole and electron-electron interactions. Interestingly enough, this 
may no longer be necessary. To illustrate this, consider the case of an s— wave pairing in a 
four-pocket model and assume for simplification that only one hole pocket is relevant to the 
pairing. Then the eigenvalue problem reduces to the set of three equations for A^, Ag, and 
Ag {Ag{k) = Ag -I- AeCOs20jfc). Solving the set, we find three solutions Xf {i — 1,2,3). In 
the absence of cos 20 terms in rij{k,p), A| = 0, and Af 2 are given by (15) with Uge — >■ '^Uee 
and — >■ 2M|g. Obviously, Uhe has to exceed a threshold, otherwise Af 2 < 0. Once the 
angle dependent terms in (9-10) become non-zero, A3 also becomes non-zero, and its sign 
depends on the interplay between oj/jg, Ogg, and /3ee- In particular, when uf^^ < UegUhh (and, 
hence, Ai_2 < 0), A3 is positive or negative depending on whether or not A > 0, where 

A = AUeeUhh {ale " /^ee) + («L + 2/3ee " Sa/igOgg) (18) 

When the angle-dependence of the electron-electron interaction can be neglected, i.e., ctge = 
Pee = 0, A3 > no matter what is the ratio of u^^ and UeeUhh- In particular, for UhhUee > u^e 
and ahe « 1, 

XI ^al >o (19) 

In other words, for one of s— wave solutions, A* > even if intra-pocket repulsions are 
the largest. The full solution of the 3 x 3 set with aee = Pee = shows that two A's 
are repulsive and one is attractive for arbitrary u\JueeUhh- When the ratio is small, the 
attractive solution is close to (19), when the ratio is large, the attractive solution is close to 
Af in (15). I illustrate this in Fig. 6 
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Nodal s± 



Modeless s± 

A. I I 

0.5 1 1.5 1 1-2 1.4 1.6 

Uu /(u u,J u, /(u u. , ) 

he ^ ee hh' he ^ ee hh' 

FIG. 6: (a) The three eigenvalues in the s— wave channel A| as functions of u'f^^/{ueeUhh) for 
Uee = Pee = and tthe = 0.4. For any u'f^^/{ueeUhh), one A| is positive (attractive), other two are 
negative. Positive A| corresponds to s± pairing. At small u\^/ {ueeUhh) pairing is induced by ahe 
and the gap has nodes on electron FSs. At large u\^/ {ueeUhh) positive Af exists already at ahe = 0, 
and the gap along electron FS has nodes only if ahe is above the threshold. The circle marks 
the area where positive and negative solutions come close to each other. The splitting between 
the two increases with ahe- (b) The regions of nodeless and nodal gap, depending on ahe and 
-u^g/MeeU/j/j. From Ref. [64]. 

There is, however, one essential difference between the cases u\^/ueeUhh > 1 and 
Uh^/ueeUhh < 1- In the first case, momentum-dependence of the interaction just modifies 
the "plus-minus" solution which already existed for momentum-independent interaction. In 
this situation, the gap along electron FS gradually acquires some cos 20 variation and re- 
mains nodeless for small ahe- In the second case, the solution with A > is induced by the 
momentum dependence of the interaction, and the eigenvalue corresponding to A3 necessary 
has Ag > Ag, i.e., s— wave gap has nodes along the electron FS [67]. In other words, the 
pairing occurs for all parameters but whether the gap is nodal or not at small ahe depends 
on the relative strength of intra-pocket and inter-pocket interactions. When intra-pocket 
interaction dominates, the gap "adjusts" and develops strong cos 20 component which does 
not couple to a momentum-independent Uee term and by this effectively reduces the strength 
of electron-electron repulsion. 

The same reasoning holds for the case of two non-equivalent hole FSs, and for 5-pocket 
models, and also for the c?— wave channel. For all cases, the solution with Aj > may exist 
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FIG. 7: Representative fits of tlie interactions rij(kp,pF) by LAHA for the 4-pocket model. Tij 
are obtained by converting the Hamiltonian, Eqs. (20), (21) from the orbital to the band basis. 
The symbols represent interactions computed numerically for the 5-band orbital model using LDA 
band structure, the black lines are the fits using Eqs. (9)-(13). The fit is for the set U = 1.67, 
J = J' = 0.21, V = 1.46, and /i = 0.08 (all in eV). A positive fi corresponds to electron doping, 
ki? in Tij(]i.p,pF) is selected along y direction on either an electron or a hole FS (its location is 
specified on top of each figure), and p_f is varied along each of FSs. The angle (p is measured 
relative to kx- 

even when intra-pocket interactions are the largest, but in this situation the gaps on the 
hole FSs have accidental nodes. The existence or non-existence of the solution at strong 
intra-pocket repulsion depends on the complex interplay between the prefactors of cos 2^ 
terms in electron-hole and electron-electron pairing vertices, see Eq. (18). 





Uhihi 




Uhih2 


Uhie 












/See 


NSF 


0.8 


0.76 


0.78 


0.46 


-0.24 


0.4 


-0.30 


0.77 


0.14 


0.09 


SF 


2.27 


2.13 


2.22 


4.65 


-0.34 


2.29 


-0.22 


3.67 


0.15 


0.04 




Uhihi 




Uhih2 


Uhie 






C^/i2e 




CKee 


Pee 


NSF 


0.7 


0.66 


-0.68 


-0.25 


-0.58 


0.24 


-0.42 


0.11 


-0.5 


0.25 


SF 


1.50 


1.40 


-1.50 


-3.73 


-0.44 


1.44 


-0.32 


1.03 


-0.49 


-0.02 



TABLE I: Table for s-wave and d— wave parameters for the same set as in Fig. 7. NSF and SF mean 
the bare interaction without the spin-fluctuation component and the full interaction, respectively. 



20 



C. How to extract rjj(k,p) from the orbital model? 



So far, in our discussion u^, a^, etc, are treated as some phenomenological inputs. 
To obtain the actual values of these parameters, one needs a microscopic model. The most 
commonly considered model for FeSCs is an effective 5-orbital model for Fe atoms with local 
intra-orbital and inter-orbital hopping integrals and intra-orbital and inter-orbital density- 
density (Coulomb) repulsions, Hund-rule exchange, and the pair hopping term. 



Hint — Uun^iMf^isi + XI ~^'^is1^it~ XI ^ J'st^^^^isa^^isa'^itaCitcr (20) 

is i,s,t^s i,s,t^s i,s,t^s f 

where riis = rii^s^ + n^^. 

The Hamiltonian Hint can be equivalently re-expressed via spin-independent interactions, 

as 

Hint ~ ^^UiiTli^s^nis^-\- ~^f^isf^it'^ ^isa^ita^^ita^ita J st ^ ^isc^isa^ita^ita 

is i,s,t^s i,s,t^s i,s,t^s f 

(21) 

where Ugt = V^t + Jstl^- 

The hopping integrals (36 total) are obtained from the fit to DFT band structure. [69] 
For the interaction parameters, the most common approximation is to assume that Ugt-, Jst 
and J'g^ are independent of the orbital indices s and as long as s 7^ t. The model can be 
also extended to include non-local Fe-Fe interactions via a pnictide [68] . 

The bare parameters in (20) and (21) are inter-related due to local spin-rotation invari- 
ance [19,35], but that invariance is broken if we view (20) and (21) as an effective low-energy 
model in which the interactions are dressed by the renormalizations coming from fermions 
with energies of order bandwidth. By this reason, in most studies t/, t/, J, and J' are 
treated as independent parameters. 

We now need to convert (20), (21) into the band basis and re-express it in the form of Eq. 
(1). This is done by transforming into the momentum space, introducing new, hybridized 
operators, which diagonalize the quadratic form, and re-expressing the interaction terms in 
(20) or (21) in terms of these new operators. The end result of this procedure is the effective 
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Hamiltonian in the band basis which has the form of Eq. (1) with rjj (k, p) given by 



where [r5f(k, k') are hnear combinations of U,U, J and J, and a.^ is the matrix element 
connecting the original fermionic operator Cp in the orbital basis with the new fermionic 
operator Oj on FS i in the band basis. The matrix elements contain information which 
orbitals mostly contribute to a particular segment of a particular FS [19,20]. Because of this, 
the interaction rjj(k, p) in the band basis generally depends on the angles along different 
FSs and contains components in all representations of the tetragonal D^^ group. 

The angle dependence of s— wave and d^^^y^ vertices agrees by symmetry with Eqs (9)- 
(13). What s a'priori unknown is how well the interactions can be approximated by the 
leading angle harmonics, i.e., whether the terms labeled as ... in (9)-(13) can actually be 
neglected. This issue was analyzed in detail in Ref . [66] , and the answer is affirmative - the 
LAHA works rather well. In Fig. 7 1 show representative fits for a particular set of parameters 
and in Table 1 1 show u^h and other parameters, extracted from the fit (in the lines marked 
marked "NSF", meaning that this is for the bare interaction, without extra spin-fluctuation 
component (see below). The results somewhat vary depending on the values of U, V, J, J', 
but in general intra-band interactions in the s— wave channel, Uee and Uhh, exceed inter band 
Uhe- This is not surprising because Uee and Uhh are essentially Coulomb interactions at small 
momentum transfers, while u^.^ is the interaction at large momentum transfer, and it should 
be smaller on general grounds. Only when V = J = J' = 0, the interaction in the band 
basis becomes independent on the momentum [70], i.e., Uee — Uhh — Uhe (this was termed 
"Coulomb avoidance" in Ref. [23]). According to Table I, intra-band interactions are also 
larger in the d-wave channel: UhihiUee > ^h-e, although the reasons why this is the case are 
not transparent. 

According to the analysis in Scc.ll, as long as intra-orbital interactions exceed inter- 
orbital interaction, Aj are negative (repulsive) if we neglect cos 2(6 terms in Fj^. If we don't 
neglect these terms, Aj can be positive even for small Uhe or Uhe but this requires that angle- 
dependent term satisfy the inequality in Eq. (18). The aeh, Oiee, and Pee, and their d— wave 
analogs in Table I are such that this condition is not me, that is for the bare interaction 
neither s— wave nor d— wave pairing is possible. 



F,,(k,k') = ^«*'*(-k)ar(k)Re[F^,^(k,k')] 



stpq 




(22) 
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The conditions on a, /3, etc depend on the parameters used for the hopings. Varying 
these parameters (e.g., making J larger) one can, in principle, met the condition in Eq. 
(18) and obtain either s— wave or d— wave solution with a positive Aj and with strong cos 2(f) 
oscillations of the gap along electron FSs. But the generic trend is that the bare interaction 
yields repulsive in s— wave and d— wave channels. 

IV. HOW TO OVERCOME INTRA-POCKET REPULSION 

How to overcome strong intra-pocket repulsion is the major issue for FeSCs. A con- 
ventional McMillan-Tolmachev renormalization [71] that reduces Coulomb repulsion and 
allows electron-phonon attraction to overcome it does not help at this stage because both 
UeeUhh and u^^ renormalize in the same way, and if repulsive UeeUhh part is initially stronger, 
the renormalization just reduces the strength of the total repulsive interaction, but doesn't 
change its sign. 

The generic idea how to overcome Coulomb repulsion goes back to Kohn-Luttinger ap- 
proach to superconductivity [32] The pairing interaction I considered so far is the bare 
interaction between a particular pair of fermions. The actual pairing interaction rjj (k, p) is 
given by the fully renormalized irreducible vertex in the particle-particle channel. This irre- 
ducible vertex is the sum of the original (bare) interaction and terms of second and higher 
order in U,V, etc. Part of these higher-order terms account for the screening of the original 
interaction by particle-hole pairs, another part can be re-cast as effective interactions medi- 
ated by collective bosonic degrees of freedom in either charge or spin channel. The issue I 
discuss now is how to obtain dressed pairing vertex. 

A. RPA approach 

One route is to adopt a semi-phenomenological approach and assume that the system is 
reasonably close to a particular density-wave ordered phase, such that the strongest compo- 
nent of the effective interaction is mediated by near-gapless fluctuations of this particular 
density- wave order. The approach goes back to Berk and Schrieffer [74], who considered 
effective pairing interaction mediated by ferromagnetic spin fluctuations. In recent years it 
has been applied to cuprates [75] and other correlated electron materials [76]. Most of 
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FIG. 8: Running vertices in SDW and SC sit channels for the 2-pocket model, panels (a)-(b), 
and the 4-pocket model, panels (c)-(d), as functions of the RG parameter L = logVF/E'. For 
the 2 pocket model A"^^^ is the largest along the whole RG flow at perfect nesting (a), but A** 
eventually wins at a finite doping when A"^^^ is cut below some E\,. For the 4-pocket model. A** 
can be the largest coupling already at perfect nesting (c), or exceed A"^^^ above some critical 
doping (d). Figures are from Ref. [64]. 

parent compounds of FeSCs possess long-range SDW order, and it is natural to assume that 
magnetically mediated pairing interaction plays the central role. 

There is no rigorously justified diagrammatic derivation of an effective fermion-fermion 
interaction mediated by collective modes, except for special cases [76]. The pragmatic 
approach is to select a set of ladder-type or bubble-type diagrams which give rise to a SDW 
order above some interaction threshold and use the same set of diagrams below the threshold. 
This approach is termed RPA by analogy with the screening problem. A generic way to apply 
RPA is to antisymmetrize the original interaction, decompose it into spin and charge channels 
and renormalize each interaction in a ladder/bubble approximation, neglecting cross-terms. 
Each 4-fermion interaction component F^^ then contains two combinations of fermionic spin 
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indices: 5a,i3Sj^s and BapB^s, where a are Pauli matrices. These two combinations have to be 
convoluted with the spin structure ia'^ of the two-particle anomalous pairing vertex. Using 
'7a'y^a,i35-t,5 = (^ps and = — Su^^ and incorporating the factor —3 into the spin part 

of the effective pairing vertex, one obtains, in the orbital basis 



where and F"^ are spin and charge components of the bare antisymmetrized interaction, 
xf^^ is the dressed static spin susceptibility and x^^^ is the dressed static charge sus- 
ceptibility. The interaction, Eq. (23), is then converted into the band basis using (22), 
and the pairing in different channels is analyzed in the same way as it was done for the 
bare interaction, i.e., either by explicitly solving integral equation on Aj in a particular 
channel [19-21,77,81], or by approximating Fjj(k, p) by the leading angular harmonics and 
solving 4 X 4 or 5 X 5 gap equations (see SecV below). The full analysis of the pairing 
problem often requires one to know the frequency dependence of but for the analysis of 
the pairing symmetry and the momentum dependence of the gap the frequency dependence 
is not overly relevant and I neglect it. 

In all model calculations, the spin susceptibility x^^^ is enhanced due to close proximity 
to SDW order, the charge susceptibility xf^"^ remains small and can be neglected. How 
(r'O^xf^^ in (23) modifies Uij depends on the type of SDW order. For moderately doped 
FeSCs X^^^{(l) is peaked at (0,7r) and (vr,0), which are the two momenta separating hole 
and electron FSs. Not surprisingly, electron-hole interaction is enhanced relative to Uee 
and Uhh- By less obvious reasons, d-wave component of electron-hole interaction is also 
enhanced, although not as strongly as Uhe- As an illustration, I show in Table I, in the lines 
marked as SF, renormalized interactions for the same set of input parameters which were 
considered in the previous section. Clearly, Uhe and Uhe are enhanced. If the enhancement 
is strong enough, cither Eq. (18) becomes valid, or uj^^/ueeUhh becomes larger than one (or 
u\^/ueeUhh becomes larger than one), and the system becomes unstable towards s-wave or 
d-wave superconductivity, whichever Aj is the largest. 

This scenario requires bare F* to be positive. Only then spin susceptibility xf^'^(g) is 
enhanced. In the orbital basis, this requirement essentially means that intra-orbital Coulomb 
repulsion U must be larger than inter-orbital U . In the band basis, the same requirement 



r?Kk,p) = [|r^ + ^r^+^(n'xr^(k-p) 




(23) 
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doping doping 



FIG. 9: Top panel: the phase diagram of the 4-pocket model at zero doping as the function of the 
strength of the angle-dependent component of the interaction, a, and Le^ = logW/ Ep, where W 
is the bandwidth and Ep is the Fermi energy. Lower panels - the phase diagrams as functions 
of doping for different Lp^. The superconducting gap in state may be nodeless or have nodes 
along electron FSs. In the region labeled as SDW, SC emerges only at a finite doping, when SDW 
instability is cut. In the region labeled as SC^, SC state emerges already at zero doping. Figure 
is from Ref. [64]. 

implies that bare Uhe must be positive [70]. What if bare Uhe is negative? Then \uhe\ niay 
still get an enhancement within RPA, but the enhancement now comes from charge (orbital) 
fluctuations [36] . The latter are strong if the system is close to a state with an orbital order. 
There are proposals [37] that broken rotational invariance observed above T/v in most of 
FeSCs [78] may be due to orbital order. [An alternative scenario [61] is an Ising order 
associated with Z2 degree of freedom separating (0, tt) and (tt, 0) spin configurations.] If 
is enhanced by orbital fluctuations such that it again exceeds UeeUhh, s-wave SC again 
develops, but the eigenvalue now corresponds to a conventional s^~^ gap, with the same sign 
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along hole and electron FSs [36]. I will compare experimental evidence for vs s"*"^ gap 
in Sec. VI. For the rest of this section I assume that > 0, i.e., the attractions is due to 
spin fluctuations. 

There are two underlying assumptions in the RPA approach. The obvious one is the 
selection of only ladder or bubble diagrams in the absence of a small parameter. The 
less obvious assumption is the very idea that the pairing can be viewed as mediated by 
collective density-wave bosonic excitations, because "mediation" implies that density-wave 
fluctuations develop at energies well above relevant energies for the pairing. This last idea is 
reasonably well justified for the cuprates, in which magnetic fluctuations develop at energies 
of a few hundred mcd, well above T,.. Whether this is also the case for FeSCs is less 
obvious as magnetism in most of these materials is itinerant (parent compounds are metals), 
and the highest Tc is at dopings at which SC and SDW ordering temperatures almost 
coincide. Whether in this situation spin fluctuations can be viewed as "pre-existing" for the 
SC problem is not obvious. It is then instructive to consider also an alternative approach 
which treats magnetic and SC fluctuations on equal footings. This alternative approach is 
based on RG technique and is termed as "KG" 

B. RG approach 

RG approach is unbiased in the sense that it doesn't assume that the pairing is mediated 
by a collecting bosonic degree of freedom. Rather, it departs from a bare Hamiltonian 
with the original 4-fermion interaction and studies how particle-particle and particle-hole 
susceptibilities evolve as one progressively integrates out contributions from fermions with 
energies larger than some running scale E. This should, in principle, address the issue 
whether spin fluctuation develop at a larger E than pairing fluctuations. In practice, this 
approach is indeed also an approximate one because even at weak coupling one cannot 
explicitly sum up contributions from fermions with energies between the bandwidth W and 
E. This can be done rigorously only if the renormalizations depend on the running scale 
logarithmically, i.e., are the functions of L = \ogW/E. Then the flow of the couphngs can 
be described by a set of differential equations dUij/dL = f{uim), where uim in f{uim) are 
renormalized interactions at the same running scale L. 

Renormalizations in a superconducting channel are indeed depend on L, but In this 
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respect FeSCs are "gifts from nature" because density-wave renormalizations involving 
fermions near hole and electron FSs also scale as logl^/£' at E above the threshold set 
by non-perfect nesting [70,79]. Then one can obtain a set of coupled RG equations for pair- 
ing and density-wave vertices. This set is called "parquet" because the renormalizations in 
the particle-particle channel are often represented as horizontal ladder diagrams, renormal- 
izations in the particle-hole (density-wave) channel as vertical ladder diagrams, and taking 
both on equal footing amounts to building perturbatively a parquet pattern. 

Prom physics perspective, the key feature of RG approach is that the pairing, SDW, 
and CDW fluctuations all "talk to each other" at intermediate energies. Not only pairing 
fluctuations are enhanced by SDW fluctuations, but also SDW fluctuations are affected by 
pairing fluctuations and so on. 

To get the idea how the RG works in FeSCs, consider momentarily a toy two-pocket model. 
Previously, in Sec. Ill A, I introduced three interactions which contribute to the pairing: 
intra-pocket hole-hole and electron-electron interactions Uhh and Ugg and inter-pocket hole- 
electron interaction Uhe > 0. Suppose that bare electron-hole interaction is weak and the 
bare A*''' < 0. Now let's add two other interactions which do not directly contribute to the 
pairing but contribute to density- wave instability: density-density and exchange interactions 
between fermions from hole and electron pockets, Udd and Uex, respectively. The coupling 
in the SDW channel is A'^^^ = Udd + Uhg. As one progressively integrates out contributions 
from high energies, all couplings evolve and obey 



Udd = 


2 1 2 
^dd + ^he 


Uex — 


'^UexiUdd - Uex) 


Uhe = 


Uhei^Udd - '^Uex 


Uhh = 


2 2 
Uhe Uf^f^ 




2 2 


Uee — 


Uhe U^e 



where the derivatives are with respect to L = (1/2) \ogW/E. We see that u^f. gets a boost 
from Udd^ which is the part of the SDW vertex, and in turn Udd gets a boost from Uhe- Solving 
the set, we find that Uhe grows under RG while Uhh and Uee decrease (and eventually change 
sign), such that below some Eq the running A*""" changes sign and becomes attractive. 

The physics behind the sign change of A*"*" is essentially the same as in RPA analysis - 
Ueh gets a boost from SDW fluctuations. But RG analysis also addresses the issue of the 
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interplay between SDW and SC orders as both are determined by the running interactions 
which "push" each other. For the two-pocket model SDW eigenvalue A'^^^ remains larger 
than A*""" along the whole RG trajectory (panels a and b in Fig. 8) , but in multi-pocket 
systems more complex behavior is possible [64,80] (see panels c and d in Fig. 8 and Fig. 
9). 

The ability to address in an unbiased way what kind of order develops first is an obvi- 
ous advantage of the RG approach. At the same same time, parquet RG has only limited 
applicability range because Eqs. (24) and the analogous equations for multi-pocket mod- 
els [64] are valid only when the running E is larger than the Fermi energy Ep- This 
wouldn't be an obstacle if the leading instability in FcSCs would develop already at E > Ep 
(Refs. [70,72,73]). But this is unlikely given that Ep ~ lO'^meV is order of magnitude larger 
than Tiv and Tc. One then has to extend the analysis to £^ < Ep. It turns out that at 
such E renormalizations in the SDW and SC channels decouple from each other [22,64] and 
each A' evolves according to d\]^/dL — (A^)^ (where now L — \ogW/E) Then SDW or SC 
instability occurs first depending on which A (either A'^-^'^ or A*""") is larger at Ep. Away 
from a perfect nesting, the logarithmical RG flow of A'^^'^ is cut below some scale Eb, while 
A** continue to grow and definitely becomes the leading instability above a certain doping. 
I illustrate this in Fig. 9. 

This consideration shows that what parquet RG actually provides are the values of the 
dressed interactions Uy at the scale of Ep. Below this scale, SDW order develops without 
input from SC channel, and on the other hand SC develops in a BCS-fashion. If the doping 
is such that SC instability is the leading instability in the problem, the symmetry and the 
structure of the SC gap is analyzed in the same way as in Sec. Ill), only Uij and Uij are 
now dressed interactions, renormalized by the RG fiow between W and Ep. From this 
perspective, RPA and RG approaches are similar in the sense that the end result of both 
approaches is the effective BCS Hamiltonian at Ep with renormalized interactions some of 
which are pushed up by SDW fluctuations. 

There is another issue with RG related to angle-dependencies of the interactions (see Sec. 
IIIC). A conventional logarithmical RG is applicable only when a and /3 terms in (9)-(13) 
are small. To order 0{a), 0{/3), angle dependencies remains intact, only the overall factors 
Uij flow [64]. Beyond this order, logarithmical and non-logarithmical terms mix up and 
the selection of diagrams can no longer be rigorously justifled. One way to proceed in this 
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FIG. 10: Representative case of small/moderate electron doping, when both hole and electron 
pockets are present. Panel a - the FS, panel b - representative fits of the interactions by LAHA 
(the dots are RPA results, the lines are LAHA expressions, Eqs (9)-(13)). Panels c and d - the 
eigenfunctions in s— wave and d— wave channels for the largest A** and A'^. From Ref. [66]. 

situation is to keep the angle dependence intact by hand and use the same equations for 
Uij as if the interactions were angle-independent [64]. Another way is to partition each FS 
into patches and write a larger set of RG equations for the interactions between different 
patches. This last approach is called functional RG (fRG) [65,80,86,87]. fRG is not an 
exact approach as it mixes logarithmical and non-logarithmical terms, but it is nevertheless 
a very powerful numerical technique to analyze the interplay between density-wave and 
superconducting instabilities in systems with angle-dependent interactions. 

Until now, the results of RPA, fRG and logarithmical RG fully agree with each other 
on (i) the interplay between SDW and SC instabilities, (ii) the symmetry of the SC state 
at different dopings, and (iii) the structure of the SC gap. This is a positive news because 
it likely implies that the underlying physics of FeSCs is quite robust. In particular, both 
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FIG. 11: The fits of the RPA interactions by LAHA and the structure of s— wave and d— wave gaps 
for the case of heavy electron doping, when only electron FSs are present. From Ref. [66]. 

fRG [80] and conventional RG [64] show that in 4-band systems SC may be the leading 
instability even without doping, while in 2-band and 5-band systems SDW is the leading 
instability at zero doping, and SC emerges only upon doping. RPA, fRG, and conventional 
RG all show that gap has nodes along electron FS in a much wider parameter range in 
4-pocket systems than in 5-pocket systems because the additional hole FS at (vr, tt) tends to 
stabilize no-nodal gap. An yet another example is the equivalence between RPA [85] and 
fRG [86] results for heavily electron doped systems, in which only electron pockets remain 
- both approaches yield a d-wave gap with no nodes. 

V. DOPING DEPENDENCE OF THE COUPLINGS 



In this section I assume that SC is the leading instability and briefly review how the 
gap symmetry and structure evolve with doping. RG, fRG, and RPA results on the doping 
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evolution of the gap agree with each other, and for definiteness I will use RPA approach 
combined with the LAHA (see Sec. III). The results differ for electron and hole dopings, 
and I consider them separately. 

A. Electron doping 

For small and moderate electron dopings, the FS consists of 4 pockets - two hole FS 
at (0,0) and two electron FSs at (0,7r) and (7r,0). Typical fits by LAHA, the parameters 
extracted from the fits, and the solutions in s-wave and d-wave channels are shown in Fig. 
10 and in Table II. It turns out [66] that some system properties are sensitive to the 
choice of the parameters, but some are quite universal. The parameter-sensitive properties 
are the presence or absence of accidental nodes in the s-wave gap (although for most of 
parameters the gap does have nodes, as in Fig. 10) and the gap symmetry itself, because for 
most of input parameters and dopings A* and A*^ remain comparable as long as both hole 
and electron FSs are present (see Table II). That d— wave state is a strong competitor in 
4-pocket systems has been first emphasized in Refs. [19,35]. The authors of [19] hinted 
that different FeSCs may have different symmetry even for the same topology of the FS. 

The universal observation is that the driving force for attraction in both s-wave and d- 
wave channels is strong inter-pocket electron-hole interaction (w/j-e and u^^e terms) no matter 
how small the hole pockets are. The gap structure actually changes only little with doping 
as long as both hole and electron pockets are present. 

TABLE II: Some of the LAHA parameters extracted from the LAHA fit in Figs. (10) and (11) for 
electron doping. Blocks (i) corresponds to Fig. (10), block (ii) corresponds to Fig. (11) (no hole 
pockets). Prom Ref [66]. 
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The situation changes qualitatively once the hole pockets disappear (Fig. 11). It is clear 
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from Table II that now the d-wave channel becomes the dominant one. Comparing the LAHA 
parameters for the two dopings, we see the reason: once the hole pockets disappear, a direct 
c?-wave electron-electron interaction becomes strong and attractive. The argument why 
this happens is as follows: [66] Uee is an antisymmetric combination of intra-pocket and inter- 
pocket electron-electron interactions Uee = ufntra ~ ""fntcr- Both Mfnter ^■^d wf^tra positive 
(repulsive), but the sign of Uee depends on the interplay between uf^^^^. and Mfntra- long as 
hole FSs are present, SF are peaked near q = (0, tt) and (tt, 0), which are an equal distance 
from the relevant momenta q = for tifn^ra ^ — i'^^ ^) '^fnter- ^^^^ situation, uf^^^.^^ 
and Mfntcr remain close in magnitude, and is small. Once the hole pockets disappear, the 
peak in the RPA spin susceptibility shifts towards (tt, tt) [85] and u^l^^j. increases more due 
to the SF component than uf^^^^. A negative wfntra ~ ""iiitcr then gives rise to a "plus-minus" 
gap on the two electron FSs. The gap changes sign under kx — >■ ky and therefore has dx2-y2 
symmetry. This pairing mechanism is essentially identical to spin-fluctuation scenario for 
d-wave pairing in the cuprates [75] . 

A different proposal has been put forward in Refs. [82] and [83]. These authors argued 
that the gap symmetry may be nodeless s— wave (equal sign of the gap on the pockets at 
(0, tt) and tt, 0). According to [82,83], s-wave pairing emerges, in some range of parameters, if 
one uses for electron-electron interaction the orbital version of the J1 — J2 model. Yet another 
proposal for strongly electron-doped FeSCs is s++ pairing driven by orbital fluctuations [84] . 



TABLE III: Some of LAHA parameters extracted from the fits in Figs. 12 and 13 for hole doping. 
Block (i) corresponds to Fig. 12 (hole and electron pockets are present), block (ii) corresponds to 
Fig. 13 ( no electron pockets). 
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FIG. 12: Representative case of small/moderate hole doping, when both hole and electron pockets 
are present. Panel a - the FS, panel b - representative fits of the interactions by LAHA (the dots 
are RPA results, the lines are LAHA expressions, Eqs (9)-(13)). Panels c and d - the eigenfunctions 
in s— wave and d— wave channels for the largest A** and A''. From Ref. [66]. 

B. Hole doping 



For small and moderate hole doping, the FS contains 5 pockets -two hole pockets at (0, 0), 
two electron pockets at (0, vr) and (vr, 0), and one more hole pocket at (vr, tt). Representative 
FSs for hole doping, typical fits by LAHA, the parameters extracted from the fit, and the 
solutions in s-wave and d-wave channels are shown in Fig. 12 and in Table III. Just like 
for electron doping, there are universal and parameter-sensitive features. The parameter- 
sensitive property is again the presence or absence of accidental nodes in the s-wave gap 
along the electron FSs, although for most of the parameters, the gap does not have nodes 
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FIG. 13: The fits of tlie RPA interactions by LAHA and the structure of s— wave and d— wave gaps 
in for strong hole doping = —0.306^), when only hole FSs are present. From Ref. [66]. 

(see Fig. 12) because the total Uhe increases once it acquires an additional contribution Uh-^e- 
There are two universal features. First, the s-wave eigenvalue is enhanced relative to a 
(i— wave one and becomes the leading instability as long as both hole and electron pockets 
are present. Second, the driving force for the attraction in both s- and d- channels is 
again strong inter-pocket electron-hole interaction {uhe and Uhe terms), no matter how small 
electron pockets are. 

The situation again changes rapidly once electron pockets disappear, see Fig. 13. Now 
electron-hole interaction becomes irrelevant, and the attractive pairing interaction may only 
be due to intra and inter-pocket interactions involving hole pockets. LAHA analysis shows 
[66] that, at least for in some range of parameters, there is an attraction in both s— wave and 
(i— wave channels, and furthermore ~ A^, see Fig. 13 The near-equivalence of s— wave and 
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d— wave eigenvalues was also found in recent unrestrictive RPA study [88]. Within LAHA, 
the attractive is due to strong intra-pocket interaction between the two hole pockets 
centered at (0,0). The s— wave gap then changes sign between these two hole pockets. 
The gap along (vr, vr) pocket is induced by a weaker inter-pocket interaction and is much 
smaller. LAHA neglects cos An(f) gap variations along the hole FSs (i.e., s— wave gaps are 
treated as angle- independent) , but the theory can indeed be extended to include these terms. 
The attractive emerges by two reasons. First, the d-wave intra-pocket interaction Uhaha 
becomes negative, second, the inter-pocket interaction Uhih^ between the two pockets at 
(0,0) becomes larger in magnitude than repulsive Uhihi and Uh2h2 (see Table III). The 
solutions with A^ > then exist separately for FSs hi^2 and for h^, the residual inter-pocket 
interaction just sets the relative magnitudes and phases between the gaps at and hi^2- 
The d-wave gap with the same structure has been obtain in the fRG analysis at large hole 
doping [87]. That study found less competition with s— wave than in PRA-based studies. 

VI. EXPERIMENTAL SITUATION 

Experimental study of FeSCs is one of the major research topics in condensed-matter 
physics since 2008 and several detailed reviews have already appeared in the literature [15, 
16,18,24,25]. Here 1 briefly discuss the experimental situation concerning the symmetry and 
the structure of the SC gap. 

As of today, there is no "smoking gun" experiment which would carry the same weight 
as phase-sensitive measurements of dx2_y2 gap symmetry in the cuprates [33]. Still, there 
is enough experimental data to minimize the number of possible gap structures. 

The theoretical proposals for the gap symmetry and structure are summarized in Sec. 
Vlll. The proposed symmetry is different for weakly/moderately doped systems with hole 
and electron FSs and for strongly doped systems where FSs of only one type are present. It 
is then instructive to consider weak/moderate and strong doping separately. 

A. Moderate doping, gap symmetry 

The candidates are s-wave (either s"^ or s++) or dx2-y2 gap. The two behave very dif- 
ferently along the hole FSs centered at (0,0) - s-wave gap is nodeless with cos 40 vari- 
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ations, while d-wave gap has nodes along kx — -^ky. ARPES measurements, both from 
synchrotron [55,89-91] and using laser hght [42], show quite convincingly that the gap 
along hole FSs is nodeless in both hole and electron-doped FeSCs. This unambiguously 
selects an s-wavc. Additional evidence in support of s-wave pairing comes from very fiat 
low-T behavior of the penetration depth in the highest Tc 1111 FeSCs systems [92]. 

B. Moderate doping, vs s++ 

The distinction between and s'^'^ gaps is a more subtle issue, particularly given that 
both belong to the same A\g representation and also because in general A\g gap on electron 
pockets may have strong oscillating component. In general, the gaps on electron and hole 
FSs have non-equal magnitudes, and the issue whether the gap is or 5++ reduces to 
whether the gap averaged over an electron FS has the same sign or opposite sign than the 
gap averaged over a hole FS. This is not a fundamental symmetry issue and, moreover, 
when cos 20 oscillations are strong, one may switch from equal to opposite signs of the 
averaged gaps by a small change of parameters [64] or by adding impurities. [93] Still, 
when oscillations are not very strong, whether the eigenfunction has or s"*"^ character is 
essential because it determines, to a large extent, whether the pairing is driven by spin or 
by orbital fluctuations (see Sec. IV). 

The experimental data most frequently cited in support of gap is the observation of 
a magnetic resonance in neutron scattering [13,94]. If, as many researchers believe, the 
resonance is a spin exciton, it exists at a momentum Q if the gaps at FS momenta kp 
and kp + Q are of opposite sign. Experimentally, the resonance is observed [13,94] near 
Q = {it, tt) in the folded BZ, which in this zone is precisely the distance between electron 
and hole FSs. The excitonic resonance then exists if the gap changes sign between hole and 
electron pockets and does not exist if the gap doesn't change sign. A similar reasoning has 
been used in identifying the the resonance seen in the cuprates with a fingerprint of d^^-y^ 
gap symmetry [95] 

From experimental perspective, the neutron peak is the resonance if it is narrow and is lo- 
cated below twice the gap value. The argument made by the supporters of s^^ scenario [36] 
is that the observed neutron peak is more broad than the resonance seen in the cuprates, 
and that there is no firm evidence that the peak energy is below 2A for the minimum gap. 
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For gap structure, there is no resonance, but there is a redistribution of the neutron 
spectral weight immediately above 2A what gives rise to a local maximum in the magnetic 
structure factor [36,96,97]. Still, the majority of researchers do believe that the observed 
neutron peak is a resonance, and the fact that it is quite broad is at least partly due to 
cos 20 gap variations along the electron FSs [97]. 

Another rather strong evidence in support of gap is the observed variation of the 
quasiparticle interference pattern in a magnetic field [98] although the interpretation of the 
data has been subject of debates [99]. It was also argued [100] that the very presence of 
the co-existence region between SC and stripe magnetism in FeSCs is a fingerprint of an 
gap, because for s~^~^ gap a first order transition between a pure magnetic and a pure SC 
state is a much more likely scenario. 

C. Moderate doping, nodal vs no- nodal gap 

Assume for definiteness that the pairing is driven by spin fluctuations and the gap has 
structure. In 2D scenario, such gap has cos 2(f) variations along electron FSs, which, according 
to theory, can be rather strong, particularly in electron-doped FeSCs. Experimental data 
show that, whether or not the gap is nodeless or has nodes, depends on the material, on the 
doping, and on whether SC co-exists with SDW order. 

1. Hole doping 

For hole-doped FeSCs (e.g. for Bai_a;Kj;Fe2As2) the data indicate that the gap is nodeless, 
away from the co-existence region. This is consistent with the theory (see Sec. V). ARPES 
experiments do not show any angular variation of the gap along both hole and electron 
FSs [55,89], but it is not entirely clear whether ARPES can at present distinguish between 
the gaps on the two electron FSs which in folded zone are both centered at (tt, tt). Thermal 
conductivity data show that k/T tends to zero in the limit of T = 0, in line with what one 
should expect for a nodeless SC [101] . Specific heat data also show non-nodal behavior [102] . 
The interpretation of the penetration depth data requires more care as the data do show a 
power-law behavior A(T) — A(0) oc T" with a ~ 2 (Refs. [103]). Such a behavior is expected 
for a SC with point nodes, but it is also expected in a wide range of T for a nodeless SC in 
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the presence of modest inter-band scattering by non-magnetic impurities [104]. Penetration 
depth measurements on artificially irradiated samples [105] support the idea that the gap 
is nodeless and power-law T" behavior of A(T) — A(0) is due to impurities. 

2. Electron doping 

For electron-doped FeSCs, e.g., 122 materials like Ba(Fei_3;Coa;)2As2 or 1111 materials 
hke NdFeAsOi_j;F-r, ARPES shows no-nodal gap along hole FS [90,91], but there are 
no data on the gap along each of the two electron FSs. At optimal doping, the data on 
both thermal conductivity [106,109] and penetration depth [109,110] are consistent with 
no-nodal gap However, the data for overdoped Ba(Fei_a;Co^)2As2 indicate that gap nodes 
may develop: the behavior of A(T) becomes more steep, and k,/T now tends to a finite 
value [106], expected for a SC with line nodes. The data also show \/H behavior of k in a 
magnetic field [106] also expected for a SC with fine nodes [107], but it was argued that 
the behavior resembhng y/H can be obtained even if gap has no nodes [108]. There 
is also clear anisotropy between in-plane conductivity and conductivity along z direction, 
what was interpreted [106] as an indication that the nodes may be located near particular 
kz- Specific heat data in overdoped Ba(Feo.9Coo.i)2As2 were also interpreted as evidence for 
the nodes. [Ill] 

The development of the nodes in gap upon electron doping is in line with the theory. 
The farther the system moves away from the SDW phase, the weaker is the increase of intra- 
band electron-hole interaction and hence the stronger is the competition from intra-band 
repulsion. As I discussed in Sec. Ill), the gap adjusts to this change by increasing its cos 20 
component in order to effectively reduce the effect of the intra-band repulsion in the gap 
equation. 

There is also experimental evidence for cos 20 gap oscillations from the observed oscilla- 
tions [112] of the field-induced component of the specific heat C{H,T) in superconducting 
FeTei_j;Se^ {x ~ 0.5). The measured C{H,T) oscillates with the direction of the applied 
field as cos 40. In theory, such an oscillation is related to the behavior of A^(0) (Ref. [62]), 
hence cos 20 gap oscillations in A lead to cos 40 oscillations in C{H, T). The observed field 
and temperature dependence of the prefactor for cos 40 term are consistent with the idea 
that the oscillations are caused by cos 20 term in A. These data were also interpreted as 
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evidence for no-nodal gap because if cos 20 gap oscillations were strong and the gap had 
nodes at accidental points, the behavior of would be more complex than the observed 
a + b cos 40. 

For LiFeAs, which is undoped but has FS structure similar to electron-doped FeSCs, 
no-nodal behavior has been observed in ARPES [113], specific heat [114], penetration 
depth [115] and NMR [43,116] measurements. The NMR data, however, were inter- 
preted [43] as evidence for a possible B-type p-wave gap {p^ + ipy). Whether this is the 
case remains an open issue. In any event, the evidence for a no-nodal gap in this material 
is overwhelming. 

3. Co-existence region 

Taken at a face value, thermal conductivity and penetration depth data indicate that 
the gap becomes nodal in the co-existence regime in both hole-doped and in electron-doped 
FeSCs. The most striking evidence comes from thermal conductivity [101,106] - in the co- 
existence regime k/T tends to a finite value at T — > and shows y/H behavior, both typical 
for a SC with line nodes. There is no theoretical understanding at present why the SC gap 
develops nodes in the co-existence regime (for a toy two-pocket model theory predicts that 
the gap should have no- nodes [117]), so this is another open issue. 

4- Isovalent doping 

Electron or hole doping is not the only way to change the properties of FeSCs. Another 
route is to replace one pnictide with the other. The most common replacement is As — >^ P. 
P-containing materials include the very first FeSC - LaFeOP, with Tc < 5K (Ref. [118]), 
the family BaFca Asi_^P^ with the highest around 30K (Ref. [119]), and LiFeP [120]. 
Penetration depth, thermal conductivity, specific heat, and NMR data [121] in these ma- 
terials all show the behavior consistent with line nodes. In particular, k scales linearly with 
T at low T and displays behavior in a magnetic field, and A(T) — A(0) is also linear in 
T down to very low T. Laser ARPES data show [42] that the gap along FS is nodeless, so 
the nodes hkely are located on electron FSs. 

On general grounds, the existence of the nodes on electron FSs is in line with theory 
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predictions particularly as BaFe2 Asi_j;Pj; has the same structure of 4 cyhndrical FSs as 
electron-doped FeSCs for which nodes are most likely. It has been argued [35] that a 
replacement of As by P changes the hight of a pnictide with respect to Fe plane, what effec- 
tively reduces intcr-pockct electron-hole interaction, in which case the gap develops nodes to 
reduce the effect of intra-pocket repulsion. However, this argument is only suggestive, and it 
is not entirely clear at the moment why all P-based FeSCs have nodes. One way to analyze 
this semi-quantitatively is to study the correlation between 2A/T'c on the hole FS and the 
presence of the nodes on electron FSs. This study shows [122] that from this perspective 
P-based FeSCs are indeed the "best case" for the gap nodes. 

Another open issue is the location of the nodes along z- direction. Oscillations of thermal 
conductivity with the direction of a magnetic field have been measured recently [123], and 
cos 40 component of these oscillations has been interpreted using the modified 2D form 
^e{kz) — Ao(l -I- a{kz)cos2(t)). The best fit to the data yields a{kz) > 1 for some and 
a{kz) < 1 for others, in which case the nodes form patches along kz- However, whether this 
is the only explanation of the data is debatable. 

I caution that although the nodes on electron FSs are most likely, they have not been 
directly observed yet, so it is still possible that the nodes are located on a hole FS, near 
particular kz-, as some of 3D theories suggest [59]. Another possibility, which is also not 
entirely ruled out, is that the system behavior near the surface, probed by ARPES, is not 
the same as in the bulk. The probability that this is the case is not high, though, because 
ARPES data are obtained using a laser light which probes states located farther from the 
surface than in conventional synchrotron-based ARPES. 

D. Strongly doped FeSCs 

1. Electron doping 

Strongly electron doped materials is recently discovered family of Aa;Fe2-j/Se2 {A — 
K,Rb,Cs) [11,12] of which K0.8Fe1.7Se2 is the most studied material. Tc in A^Fe2_ySe2 
is rather high, almost 40K. ARPES shows [12] that only electron FSs are present in 
A3;Fe2-ySe2, while hole pockets are at least 60mey from the FS, although hole dispersion 

above GOmeV is still clearly visible in ARPES. Two electron FSs are at (0,7r) and (tt, 0), 
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like in other FeSCs, and there is, possibly, another electron FS at (0,0). RPA, LAHA and 
fRG calculations for these systems predict that the gap should have a d-wave symmetry, at 
least for the case when the FSs are only at (0,7r) and (vr, 0). A d-wave symmetry in this 
situation means that the gaps on the two electron FSs behave as Ao(±l + a cos 20), and all 
calculations yield a < 1, i.e., no nodes (neglecting 3D effects). The theoretical alternative 
is s"*""*" symmetry by one reason [84] or the other [82,83]. At present, both ARPES and 
specific heat data point that the gap is nodeless, at least for most of values, but whether 
the gap is an s— wave or a nodeless d— wave remains to be seen. 

2. Hole doping 

Strongly hole doped system is KFe2As2 (Tc = 3K) which is at the opposite end from 
parent BaFe2As2 in the family of K2.Bai_a;Fe2As2. According to ARPES [40], this system 
has no electron pockets and has hole pockets at (0,0). There may also be additional hole 
pockets, but this is not entirely clear. Both thermal conductivity and penetration depth 
measurements clearly point to nodal behavior [41]. There is, however, no "smoking gun" 
symmetry-sensitive measurement, so whether the gap is a d-wave or an s-wave with nodes 
due to strong cos 40 gap component on one of the FSs remains an open issue 

E. Summary 

Overall, the agreement between theory and experiment with respect to gap symmetry and 
structure is reasonably good. Theory predicts that the gap symmetry in weakly and moder- 
ately doped FeSCs is an s-wave. This is consistent with ARPES data. Many theorist argue 
that an s-wave gap changes sign between hole and electron FSs. Quasiparticle interference 
and neutron scattering data are consistent with this picture, if, indeed, the neutron peak is 
the resonance. Further, theory predicts that gap has oscillations along electron FS, and 
these oscillations give rise to accidental nodes, which are more likely for systems with 2 hole 
and 2 electron cylindrical FSs than in systems with an additional hole cyhndrical FS. This 
is also generally consistent with the experiments, although at this moment the agreement 
is more qualitative than quantitative. For strongly electron-doped FeSCs, the theory based 
on RPA, LAHA, and fRG predicts a nodeless d-wave gap, except, perhaps, the region near 
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kz — 7r/2. This is neither confirmed not disproved by the experiments. For strongly hole- 
doped FeSCs, with only hole pockets, experiments point to the presence of the nodes, but 
whether the gap is a c?— wave or a nodal s— wave has not been firmly established. Finally, in 
the co-existence region of SC + SDW, theory prediction for a 2-band toy model is a nodeless 
gap. This is consistent with the data for weak SDW order, but apparently inconsistent 
with the data for FeSCs with strong SDW order, which show behavior consistent with the 
nodes in the SC gap. What is the gap structure deep in the co-existence region remains to 
be seen. 



VII. CONCLUSION 



The analysis of the gap symmetry and structure in FeSCs is a fascinating subject because 
of multi-orbital/multi-band nature of these materials. It is now well understood that in 
multi-band systems a conventional notion that s-wave gap is nodeless, d-wave has 4 nodes, 
etc, does not work, and s-wave gap may have nodes, while d-wave gap may remain nodeless. 
Furthermore, nodes may appear or disappear, depending on doping and other external 
conditions. Another peculiarity of FeSCs is a close proximity between pairing eigenvalues Aj 
from different symmetry representations, e.g., Aig and Big. Because of proximity, even the 
gap symmetry may change upon doping or other external perturbation. Also, below T^, the 
system may lower the energy by moving into a mixed state with, e.g., s + id gap function. 
To analyze this issue, one needs to go beyond the analysis presented in this review and solve 
the full non-linear gap equation. 

The existence of superconductivity at the "end points" of the phase diagram, when only 
hole pockets or only electron pockets are present, is another fascinating issue. There is 
reasonably good understanding of a d-wavc pairing at strong hole doping, but there are 
competing proposals for the gap symmetry at strong electron doping, and so far experiments 
cannot distinguish between these proposals. Finally, the gap structure in the co-existence 
regime with SDW is yet another fascinating unexplored issue, and comprehensive analysis 
of SC in the co-existence phase is clearly called for. 
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VIII. SUMMARY POINTS 

• For weakly and moderately electron-doped FeSCs, and dx2_y2 pairing channels are 
nearly degenerate. The driving force for the pairing in both channels is inter-pocket 
electron-hole interaction, enhanced by spin fluctuations, d-wave gap has nodes at 
kx — ±ky on the hole FSs, gap has cos 2(f) variations on electron FSs and may have 
nodes. The probability that gap has nodes increases with electron doping. 

• For strongly electron doped FeSCs, when only electron pockets remain, some results 
show that the gap has a d— wave symmetry, which in this situation implies that angle- 
independent component of the gap changes sign between the two electron pockets. 
The driving mechanism for the pairing is a direct d— wave attraction between electron 
pockets, again enhanced by spin fluctuations. The gap is nodeless, but may acquire 
nodes near kz = 7r/2 when 3D effects and/or hybridization between electron pockets 
are included. A competing theory proposal, based on orbital J — 1 — J — 2 model, is 
that the gap is sign-preserving s-wave. 

• For weakly and moderately hole-doped FeSCs, s"^ pairing is the dominantinstability. 
The driving force is again inter-pocket electron-hole interaction, enhanced by spin 
fluctuations. The gap has cos 20 variation along electron FSs but hkely does not have 
nodes. 

• For strongly hole-doped FeSCs, when only hole FSs remain, the gap symmetry may 
be d— wave, with nodes at kx — -^ky on all hole FSs. The d-wave pairing mechanism is 
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the combination of a d-wave attraction within the hole pocket at (tt, tt) and a strong 
intra-pocket d-wave interaction between the two hole pockets centered at (0,0). A 
competing theory proposal is that the gap is s-wave, with potential nodes due to 
strong cos 40 variation along hole FSs centered at (0, 0). 

• For real 3D systems, one scenario is that the 2D picture survives, and nodes, if present, 
are "vertical" nodes, located on electronic FSs at given and ky but arbitrary k^. 
Another scenario is that nodes on electron FSs exists only in some ranges of kz- And 
the third, more radical scenario, is that the nodes appear as "horizontal" nodes near 
particular kz, either on electron or on hole FSs. 

• An alternative scenario for all FeSCs is that the gap is a conventional, sign-preserving 
s"*""*" gap, but still may have cos 20 oscillations on electron FSs. Such a gap appears if 
inter-pocket electron-hole interaction Uhe is negative (attractive) and is enhanced by 
charge (orbital) fluctuations and by phonons 

IX. A LIST OF FUTURE ISSUES 

• What would be the "smoking gun" experiment to distinguish between and s^^ 
gap? 

• Why all P-based FeSCs appear to have nodes, and where these nodes are located? 

• What is the structure of the SC gap in the co-existence region with SDW order? 

• What is the pairing mechanism and the gap symmetry and structure in heavily 
electron-doped Kj;Fe2_ySe2 and in heavily hole-doped KFe2As2? 

• Is there a mixed s + id SC state in FeSCs at T = 0? 

• Is there a FeSC with a p-wave symmetry? 
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